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Extensive systematizations of theoretical and experimental nuclear densities and of optical po- 
tential strengths extracted from heavy-ion elastic scattering data analyses at low and intermediate 
energies are presented. The energy-dependence of the nuclear potential is accounted for within a 
model based on the nonlocal nature of the interaction. The systematics indicates that the heavy-ion 
nuclear potential can be described in a simple global way through a double- folding shape, which 
basically depends only on the density of nucleons of the partners in the collision. The possibility of 
extracting information about the nucleon-nucleon interaction from the heavy-ion potential is inves- 
tigated. 
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1. Introduction 

The optical potential plays a central role in the description of heavy-ion collisions, since it is widely used in studies 
of the elastic scattering process as well as in more complicated reactions through the DWBA or coupled-channel 
formalisms. This complex and energy-dependent potential is composed of the bare and polarization potentials, the 
latter containing the contribution arising from nonelastic couplings. In principle, the bare (or nuclear) potential 
between two heavy ions can be associated with the fundamental nucleon-nucleon interaction folded into a product 
of the nucleon densities of the nuclei JjJ. Apart from some structure effects, the shape of the nuclear density along 
the table of stable nuclides is nearly a Fermi distribution, with diffuseness approximately constant and radius given 
roughly by R = tq A 1 / 3 , where A is the number of nucleons of the nucleus. Therefore, one could expect a simple 
dependence of the heavy- ion nuclear potential on the number of nucleons of the partners in the collision. In fact, 
analytical formulae have been deduced [|]-|| for the folding potential, and simple expressions have been obtained at 
the surface region. An universal (system-independent) shape for the heavy-ion nuclear potential has been derived 0] 
also in the framework of the liquid-drop model, from the Proximity Theorem which relates the force between two 
nuclei to the interaction between flat surfaces made of semi-infinite nuclear matter. The theorem leads Q to an 
expression for the potential in the form of a product of a geometrical factor by a function of the separation between 
the surfaces of the nuclei. 

The elastic scattering is the simplest process that occurs in a heavy-ion collision because it involves very litle 
rearrangement of matter and energy. Therefore, this process has been studied in a large number of experimental 
investigations, and a huge body of elastic cross section data is currently available. The angular distribution for elastic 
scattering provides unambiguous determination of the real part of the optical potential only in a region around a 
particular distance Q hereafter referred as the sensitivity radius (i?s)- At energies close to the Coulomb barrier the 
sensitivity radius is situated in the surface region. In this energy region, the systematization J?]]^] of experimental 
results for potential strengths at the sensitivity radii has provided an universal exponential shape for the heavy-ion 
nuclear potential at the surface, as theoretically expected, but with a diffuseness value smaller than that originally 
proposed in the proximity potential. 

In a recent review article || the phenomenon of rainbow scattering was discussed, and it was emphasized that the 
real part of the optical potential can be unambiguously extracted also at very short distances from heavy-ion elastic 
scattering data at intermediate energies. Such a kind of data has been first obtained for a-particle scattering from a 
variety of nuclei over a large range of energies |j|-[n|, and later for several heavy- ion systems. However, differently from 
the case for the surface region (low energy), a systematization of potential strengths at the inner distances has not 
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been performed up to now, probably because the resulting phenomenological interactions have presented significant 
dependence on the bombarding energies. Several theoretical models have been developed to account for this energy- 
dependence through realistic mean field potentials. Most of them are improvements of the original double-folding 
potential with the nucleon-nucleon interaction assumed to be energy- and density-dependent || . Another recent and 
successful model [^2|-[l4| associates the energy-dependence of the heavy-ion bare potential with nonlocal quantum 
effects related to the exchange of nucleons between target and projectile, resulting in a very simple expression for the 
energy-dependence of the nuclear potential. Using the model of Refs. [|12|-[l4|, in the present work we have realized 
a systematization of potential strengths extracted from elastic scattering data analyses, considering both: low (near- 
barrier) and intermediate energies. The systematics indicates that the heavy-ion nuclear potential can be described 
in a simple global way through a double-folding shape, which basically depends only on the number of nucleons of 
the nuclei. 

The paper is organized as follows. In Section 2, as a preparatory step for the potential systematization, an extensive 
and systematic study of nuclear densities is presented. This study is based on charge distributions extracted from 
electron scattering experiments 151^| as well as on theoretical densities derived from the Dirac-Hartree-Bogoliubov 
model |Ti| . In Section 3, analytical expressions for the double-folding potential are derived for the whole (surface 
and inner) interaction region, and a survey of the main characteristics of this potential is presented. Section 4 
contains the nonlocal model for the heavy-ion bare interaction, including several details that have not been published 
before. Section 5 is devoted to the nuclear potential systematics. In Section 6, we discuss the role played by the 
nucleon-nucleon interaction, and we present, in a somewhat speculative way, an alternative form for the effective 
nucleon-nucleon interaction, which is consistent with our results for the heavy-ion nuclear potential. Finally, Section 
7 contains a brief summary and the main conclusions. 



2. Systematization of the nuclear densities 

According to the double-folding model, the heavy-ion nuclear potential depends on the nuclear densities of the nuclei 
in collision. Thus, a systematization of the potential requires a previous systematization of the nuclear densities. In 
this work, with the aim of describing the proton, neutron, nucleon (proton+neutron), charge and matter densities, we 
adopt the two-parameter Fermi (2pF) distribution, which has also been commonly used for charge densities extracted 
from electron scattering experiments [|l5|| . The shape, Eq. 1 and Fig. 1, of this distribution is particularly appealing 
for the density description, due to the flatness of the inner region, that is associated with the saturation of the nuclear 
medium, and to the rapid fall-off (related to the diffuseness parameter a) that brings out the notion of the radius, 
i?o> of the nucleus. 

P ^ = 7T P °{r-R \ (!) 

1 + exp (J) 

The a and ^0 parameters are connected by the normalization condition: 

p DO 

4tt / p(r) r 2 dr = X , (2) 







where X could be the number of protons Z, neutrons N or nucleons A = N + Z. In our theoretical calculations, the 
charge distribution (p c h) has been obtained by folding the proton distribution of the nucleus (p p ) with the intrinsic 
charge distribution of the proton in free space (p c hp) 



Pch(r) = J p p (r') Pch P (r- r') dr 1 , (3) 

where p c h P is an exponential with diffuseness a c h p = 0.235 fm. In an analogous way, we have defined the matter 
density of the nucleus by folding the nucleon distribution of the nucleus with the intrinsic matter distribution of the 
nucleon, which is assumed to have the same shape of the intrinsic charge distribution of the proton. For convenience, 
the charge and matter distributions are normalized to the number of protons and nucleons, respectively. 

In order to systematize the heavy- ion nuclear densities, we have calculated theoretical distributions for a large 
number of nuclei using the Dirac-Hartree-Bogoliubov (DHB) model p7| . The DHB calculations were performed using 
the NL3 parameter set fis|| . This set was obtained by adjusting the masses, and the charge and neutron radii of 10 
nuclei in the region of the valley of stability, ranging from 16 to 214 Pb, using the Dirac-Hartree-BCS (DH-BCS) 
model. For the cases in which they have been performed, calculations using this parameter set and either the DHB |17| ] 
or the DH-BCS |lq-G(J] model have shown very good agreement with experimental masses and radii. In the present 



work, we have also used the results of previous systematics for charge distributions (T^JT^], extracted from electron 
scattering experiments, as a further check of our DHB results. All the theoretical and most of the "experimental" 
densities are not exact Fermi distributions. Thus, with the aim of studying the equivalent diffuseness of the densities, 
we have calculated the corresponding logarithmic derivatives (Eq. 4) at the surface region (at r ~ i?o + 2 fm). 

a ~-^ ( 4 ) 
dr 

Fig. 2a shows the results for the experimental charge distributions: the diffuseness values spread around an average 
diffuseness a c = 0.53 fm, with standard deviation 0.04 fm. Most of this dispersion arises from experimental errors. 
Indeed, we have verified that different analyses (different electron scattering data set or different models for the charge 
density) for a given nucleus provide diffuseness values that differ from each other by about 0.03 fm. Therefore, the 
experimental charge distributions are compatible, within the experimental precision, with a constant diffuseness value. 
The theoretical charge distributions present similar behavior (Fig. 2b), with average value slightly smaller than the 
experimental one. In this case, the observed standard deviation, 0.02 fm, is associated with the effects of the 
structure of the nuclei. Despite the trend presented by the neutron and proton diffuseness (Fig. 2c), all the nucleon 
distributions result in very similar diffuseness values (a_/v = 0.48 /to), with standard deviation 0.025 fm. Due to the 
folding procedure, the matter distributions present diffuseness values significantly greater (om = 0.54 /to) than those 
for the nucleon distributions. Taking into account that the theoretical calculations have slightly underestimated the 
experimental charge diffuseness, we consider that more realistic average values for the nucleon and matter density 
diffuseness are 0.50 fm and 0.56 /to, respectively. A dispersion (<r ) of about 0.025 fm around these average values 
is expected due to effects of the structure of the nuclei. 

The root-mean-square (RMS) radius of a distribution is defined by Eq. 5: 



/ J> 2 p{r) dr 

r ^-\j j p{r)d f ■ W 

We have determined the radii Rq for the 2pF distributions assuming that the corresponding RMS radii should be 
equal to those of the experimental (electron scattering) and theoretical (DHB) densities. The results for Ro from 
theoretical charge distributions (Fig. 3b) are very similar to those from electron experiments (Fig. 3a). This fact 
indicates that the radii obtained through the theoretical DHB calculations are quite realistic. The nucleon and matter 
densities give very similar radii (Fig. 3d), which are well described by the following linear fit: 

Ro = 1.31 A 1/3 - 0.84 fm. (6) 

Due to effects of the structure of the nuclei, the Rq values spread around this linear fit with dispersion ctr = 0.07 fm, 
but the heavier the nucleus is, the smaller is the deviation. In Fig. 4 are shown the theoretical (DHB) nucleon densities 
for a few nuclei, and the corresponding 2pF distributions with a = 0.50 fm and Rq values obtained from Eq. 6. 



3. Essential features of the folding potential 

The double-folding potential has the form 

V F (R) = J Pi{n) P2{r 2 ) v NN (R -fi+fz) df{ df 2 , (7) 

where R is the distance between the centers of the nuclei, pi are the respective nucleon distributions, and VNN{r) is 
the effective nucleon-nucleon interaction. The success of the folding model can only be judged meaningfully if the 
effective nucleon-nucleon interaction employed is truly realistic. The most widely used realistic interaction is known 
as M3Y jjjfjl, which can usually assume two versions: Reid and Paris. 

For the purpose of illustrating the effects of density variations on the folding potential, we show in Fig. 5 the 
results obtained for different sets of 2pF distributions. In Section 2, we have estimated the dispersions of the Ro and 
a parameters, o\r w 0.07 and a a w 0.025 fm, that arise from effects of the structure of the nuclei. Observe that these 
standard deviations are one half of the corresponding variations considered in the example of Fig. 5, ARq = 0.14 fm 
and Aa = 0.05 fm. The surface region of the potential (R > R\ + R 2 ) is much more sensitive to small changes of the 
density parameters than the inner region. Our calculations indicate that, due to such structure effects, the strength 
of the nuclear potential in the region near the barrier radius may vary by about 20%, and the major part of this 



variation is connected to the standard deviation of the parameter a. Therefore, concerning the nuclear potential, the 
effects of the structure of the nuclei are mostly present at the surface and mainly related to the diffuscness parameter. 

The six-dimensional integral (Eq. 7) can easily be solved by reducing it to a product of three one-dimensional 
Fourier transforms Q, but the results may only be obtained through numerical calculations. In order to provide 
analytical expressions for the folding potential, we consider, as an approximation, that the range of the effective 
nucleon-nucleon interaction is negligible in comparison with the diffuseness of the nuclear densities. In this zero-range 
approach, the double-folding potential can be obtained from: 



v NN (r) » V 5(r) V F {R) = 
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As discussed in Section 2, the heavy-ion densities involved in Eq. 8 are approximately 2pF distributions, with Rq ^> a. 
In the limit a — > 0, the double-integral results in 



V F (R<R 2 -R 1 ) = V poi P02 gTri?? , 



(9) 



V F {R 2 - Ri < R < Ri + R 2 ) = V poi p 02 - 
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(10) 



V F (R >Ri + R 2 ) = 0, 



(11) 



where s = R — + R 2 ), 1Z — 2R 1 R 2 /(Ri + R 2 ), ( — 1Z/{R\ + R 2 ), r = s/lZ, R\ and R 2 are the radii of the nuclei 
(hereafter we consider R 2 > R\). We need a further approximation to obtain analytical expressions for the folding 
potential in the case of finite diffuseness value. 

The Fermi distribution may be represented, with precision better than 3% for any r value (see Fig. 1), by: 



r-Rr, 



(12) 



C (x<0) = l-~e* + ~e 2 * 



(13) 



1 - -e 



(14) 



This approximation is particularly useful in obtaining analytical expressions for integrals that involve the 2pF distri- 
bution. If both nuclei have the same diffuseness a, the double-integral (Eq. 8) can be solved analytically using the 
approximation represented by Eq. 12, and the result expressed as a sum of a large number of terms, most of them 
negligible for a <C Rq. Rather simple expressions can be found after an elaborate algebraic manipulation: 
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V F (R >Ri + R 2 ) w V poi P02 ™ 2 K g(r) f{s/a) , 



(17) 



with rj = a/7?., e = s/a. The functions g and / are given by: 



f(s/a) = (1 + s/a) e- s ' a 



(19) 



If the nuclei have slightly (about 10%) different diffuscncss, the formulae are still valid with a « {a\ + 02) /2. As an 
example of the precision of the analytical expressions above, we exhibit in Fig. 6 the results of numerical calculations 
(Eq. 8) and compare them with those from Eqs. 15, 16 and 17, and also with the exact expressions for a — 0, Eqs. 
9, 10 and 11. 

Eq. 17 presents some similarity with the proximity potential g): 

V P = 2-KTTZd $(s/d) , (20) 

where d is the "surface width" , and $ is an universal (system-independent) function. For a 2pF distribution, the 
surface width is related to the diffuseness parameter through: d ~ (n/ a p3|. The theoretical value adopted for d 
is 1 fm |J], which corresponds to a diffuseness a rs 0.55 fm. Taking into account that the T value is rather system- 
independent [El , systematizations of heavy-ion potential strengths extracted from elastic scattering data analyses have 
been performed by using the following expression, which should be valid for surface distances, 

= V, e-'° . (21) 

The resulting experimental a values are quite similar, a ~ 0.62 fm but smaller than the theoretical prediction 
of the proximity potential a ~ 0.75 fm |5|. Such systematics have included only experimental potential strengths in 
the surface region, in contrast to the case of the proximity potential where V/1Z should be an universal function of s 
also for inner distances. The proximity potential is not fully agreeing with our results for the double-folding potential 
in the zero-range approach (see Fig. 7). In fact, Eq. 17 indicates that a better choice for an universal quantity at the 
surface region would be 

Vred(s > 0) = VP , (22) 

P01 P02 TTa^Tl g[T) 

which results (from Eqs. 17, 19 and 22) the system-independent expression 

V red {s> 0)^V (1 + s/a) e- s / a . (23) 

However, it is not clear that one can find a simple form for such a universal quantity at inner distances from Eqs. 15 
and 16. In Section 5, the reduced potential, V re d, is useful for addressig the potential strength systematization. Thus 
we define V re d for s < through the following trivial form: 

V red (s < 0) = V . (24) 

The end of this section is devoted to the study of the effect on the folding potential of a finite range for the effective 
nucleon-nucleon interaction. The tri-dimensional delta function, Vq S(f), can be represented through the limit a — > 
applied to the finite-range Yukawa function 

YAr)=V ^^. (25) 

Fig. 8 shows a comparison of folding potentials in the zero-range approach (Eq. 8) with the result obtained (from 
Eq. 7) using an Yukawa function for the effective nucleon-nucleon interaction. The finite range is not truly significant 
at small distances, and can be accurately simulated at the surface, within the zero-range approach, just by slightly 
increasing the diffuseness of the nuclear densities. 



4. A nonlocal description of the nucleus-nucleus interaction 

Before proceeding with the systematization of the potential, we first set the stage for the model of the heavy-ion 
nuclear interaction |l2]-|l4|]. When dealing with nonlocal interactions, one is required to solve the following integro- 
differential equation 



- —V 2 f(R) + [V C {R) + V pol {R,E) +iW pol (R,E)}V{R) + J U(R,R') ^{R') dR' = E^(R) . (26) 

Vc is the Coulomb interaction assumed to be local. V po i and W po i are the real and imaginary parts of the polarization 
potential, and contain the contribution arising from nonelastic channel couplings. The corresponding nonlocality, 
called the Feshbach nonlocality, is implicit through the energy-dependence of these terms, consistent with the disper- 
sion relation [p2| . U(R, R') is the bare interaction and the nonlocality here, the Pauli nonlocality, is solely due to the 
Pauli exclusion principle and involves the exchange of nucleons between target and projectile. 

Guided by the microscopic treatment of the nucleon- nucleus scattering ]2^-|27|| , the following ansatz is assumed for 
the heavy-ion bare interaction |l3[ ] 



U{R, R') = Vi 
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(27) 



where b is the range of the Pauli nonlocality. Introduced in this way, the nonlocality is a correction to the local model, 
and in the b — > limit Eq. 26 reduces to the usual Schroedinger differential equation. The range of the nonlocality 
can be found through b w bomo/n PSfl i where bo = 0.85 fm is the nucleon-nucleus nonlocality parameter J23|| , mo 
is the nucleon mass, and (i is the reduced mass of the nucleus-nucleus system. This type of very mild nonlocality in 
the nucleon-nucleus and nucleus-nucleus interaction is to be contrasted with the very strong nonlocality found in the 
pion-nucleus interaction in the A-region p9[ ]. In such cases, even the concept of an optical potential becomes dubious. 
In our case, however, we are on very safe ground. 

The relation between the nonlocal interaction and the folding potential is obtained from |H| 

V NL (R) = V F {R) . (28) 

Due to the central nature of the interaction, it is convenient to write down the usual expansion in partial waves, 
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where Qe are polynomials and (f> is the angle between R and R' 23 1 . Thus, the integro-differential equation can be 
recast into the following form 



h z d 2 u e (R) 
27 dR 2 
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V e (R, R') ue(R') dR' . (32) 



When confronting theory and experiment, one usually relics on the optical model with a local potential. This brings 
into light the issue of extracting from Eq. 32 a local- equivalent (LE) potential 



Vle(R,E) + iWle{R, E) 
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u e (R) Jo 
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(33) 



The presence of the wave-function in Eq. 33 indicates that the LE potential is complex and also t- and energy- 
dependent. Despite its complex nature, the LE potential is not absorptive, ( v I'|T / 1 / l,e| , I') = 0; this statement can 
be demonstrated by considering that the nonlocal interaction is real and symmetrical, Ve(R, R') = Vg(R',R). For 
neutron- nucleus systems, the LE potential is only weakly ^-dependent, and an approximate relation to describe its 
energy-dependence has been obtained . A generalization of this relation for the ion-ion case is given by (lj,[l]| : 



V LE (R,E)^V F (R) e 



- 1 [E-Vc{R)-V LE (R,E)} 



(34) 



with 7 = /ib 2 /2h 2 . In order to provide an example of the precision of expression 34, in Fig. 9 the corresponding 
result is compared to the exact LE potential (Eq. 33) obtained from the numerical resolution [Ol of the respective 



integro-differential equations (Eq. 32). The local-equivalent potential is quite well described by Eq. 34 for any I 
value, except at very small distances (R rs 0) that are not probed by heavy-ion experiments. 

Expression 34 has accounted for the energy-dependence of experimentally extracted potential strengths for several 
systems in a very large energy range ]l2|-^4|. At near-barrier energies, E sa Vc(Rb) + Vle(Rb), the effect of the 
Pauli nonlocality is negligible and Vle(R,E) ss Vf(R), but the higher the energy is, the greater is the effect. At 
energies about 200 MeV/ nucleoli the local-equivalent potential is about 1 order of magnitude less intense than the 



corresponding folding potential (see examples in Refs. |12|l3j). In a classical physics framework, the exponent in Eq 



34 is related to the kinetic energy (Ek) and to the local relative speed between the nuclei (v) by 

v 2 = -E k (R) = - [E-Vc(R) - V LE (R,E)} ■ (35) 

and Eq. 34 may be rewritten in the following form 

V LE (R,E) « V F (R) e -[^ob v/(2h)f M Vp ^ e -4^/c 2 ; ( 36 ) 

where c is the speed of light. Therefore, in this context the effect of the Pauli nonlocality is equivalent to a velocity- 
dependent nuclear interaction (Eq. 36). Another possible interpretation is that the local-equivalent potential may 
be associated directly with the folding potential (Eq. 37), with an effective nucleon-nucleon interaction (Eq. 38) 
dependent on the relative speed (v) between the nucleons 



V LE {R,E) =Vf = J Pi(n) P2{r2) v NN (v,R-fi +f 2 ) dfx df 2 , (37) 
v NN { Vl f)=v f {r)e-^' c2 . (38) 

5. The systematization of the nuclear potential 

As already mentioned, the angular distribution for elastic scattering provides an unambiguous determination of the 
real part of the optical potential in a region around the sensitivity radius {Rs)- For bombarding energies above (and 
near) the barrier, the sensitivity radius is rather energy-independent and close to the barrier radius (Rb), while at 
intermediate energies much inner distances are probed. At sub-barrier energies, the Rs is strongly energy-dependent, 
with its variation connected to the classical turning point; this fact has allowed the determination of the potential in 
a wide range of near-barrier distances, Rb < Rs < Rb + 2 fm. With the aim of avoiding ambiguities in the potential 
systematization, we have selected "experimental" (extracted from elastic scattering data analyses) potential strengths 
at the corresponding sensitivity radii, from works in which the Rs has been determined or at least estimated. In 
several articles, the authors claim that their data analyses at intermediate energies have unambiguously determined 
the nuclear potential in a quite extensive region of interaction distances. In such cases, we have considered potential 
strength "data" in steps of 1 fm over the whole probed region. Tables I and 2 provide the systems included in the 
nuclear potential systematics for the sub-barrier and intermediate energies, respectively. For the energy region above 
(and near) the barrier, the present systematics contains potential strengths for a large number of different heavy-ion 
systems from the previous Christensen and Winther's systematization Q. Our systematics is not even near to being 
complete, but it is rather extensive and diversified enough to account well for the very large number of data that have 
been obtained in the last decades. 

The experimental potential strengths represent the real part of the optical potential, which corresponds to the 
addition of the bare and polarization potentials. The contribution of the polarization to the optical potential depends 
on the particular features of the reaction channels involved in the collision, and is therefore quite system-dependent. 
If this contribution were very significant, it would be too difficult for one to set a global description of the heavy- 
ion nuclear interaction. In the present work, we neglect the real part of the polarization potential and associate the 
experimental potential strengths (V E xp) with the bare interaction (Vle)- The success of our findings seems to support 
such a hypothesis. 

In analysing experimental potential results for such a wide energy range and large number of different systems, 
we consider quite appropriate the use of system- and energy-independent quantities. We have removed the energy- 
dependence from the experimental potential strengths through the calculation of the corresponding folding potential 
strengths, V F ~Exp, based on Eq. 34. The system-dependence of the potential data set has then been removed with 
the use of the experimental reduced potential, V re d-Exp- For s > this quantity was calculated from Eq. 22, and for 
inner s values we have adopted the following simple definition 



V re d-Exp = Vo TT — , (39) 

VF-Teo 

with VF-Teo calculated through Eq. 8. The other useful quantity is the distance between surfaces: s = Rs— (R1+R2), 
where Rs is associated to the sensitivity radius, and the radii of the nuclei have been obtained from Eq. 6. 

In Fig. 10 (top), the experimental reduced potential strengths are confronted with the theoretical prediction (Eqs. 
23 and 24) for different diffuseness values. The fit to the data in the inner region (s < 0) results unambiguously in the 
value Vo = —456 MeV fm , and is quite insensitive to the diffuseness parameter, in agreement with the discussion 
about the folding features of Section 3. The fit for s > is sensitive to both: Vb and a, and the corresponding best fit 
values are a = 0.56 fm and the same Vq found for the inner region. The standard deviation of the data set around 
the best fit (solid line in Fig. 10 - top) is 25%, a value somewhat greater than the dispersion (20%) expected to 
arise from effects of the structure of the nuclei (as discussed in Section 3). We believe that the remaining difference 
comes from two sources: uncertainties of the experimentally extracted potential strengths and the contribution of 
the polarization potential that we have neglected in our analysis. We point out that the best fit diffuseness value, 
a = 0.56 fm, is equal to the average diffuseness found (Section 2) for the matter distributions and greater than the 
average value (a = 0.50 fm) of the nucleon distributions. This is a consistent result because we have calculated 
the reduced potential strengths based on the zero-range approach (through Eqs. 8, 22, 23 and 24). As discussed 
in Section 3, the effect of a finite-range for the effective nucleon-nucleon interaction can be simulated, within the 
zero-range approach, by increasing the diffuseness of the (nucleon) densities of the nuclei. This subject is dealt with 
more deeply in the next Section. 

In order to characterize the importance of the Pauli nonlocality, in Fig. 10 (bottom) are shown the results for the 
reduced potential through calculations performed without the correction (Eq. 34) due to the energy-dependence of 
the LE potential, i.e. associating the experimental potential strengths directly with the folding potential. The quality 
of the corresponding fit (Fig. 10 - bottom) is similar to that obtained with the nonlocality (Fig. 10 - top), but the Vo 
and a parameters are significantly different. In the next Section, we show that the values found without considering 
the nonlocality, a = 0.61 fm and Vq — —274 MeV fm 3 , seem to result in an unrealistic nucleon-nucleon interaction. 



6. The effective nucleon-nucleon interaction 

After removing the energy-dependence of the experimental potential strengths, the corresponding results are com- 
patible with the double- folding potential in the zero-range approach (Eq. 8), provided that the matter densities 
of the nuclei be adopted in the folding procedure instead of the nucleon densities. In this section, we study the 
consistency of our results for the nuclear potential in the case that the double-folding model is treated in the more 
common interpretation: the nucleon distributions and a finite-range nucleon-nucleon interaction are assumed in Eq. 

7. With the purpose of keeping the comparison between experimental and theoretical results through the use of 
system-independent quantities, it is necessary to change the definition of the experimental reduced potential 

V re d-Exp = V re d-Teo ^~ ^ , (40) 

VF-Teo 

where Vp-Teo is now calculated through Eq. 7. V re d-Teo is still obtained from Eqs. 23 and 24, with the Vo parameter 
being associated to the volume integral of the effective nucleon-nucleon interaction (actually, this same procedure has 
also been adopted in the zero-range case) 



V = 47T / v NN (r) r z dr . (41) 

The effective nucleon-nucleon interaction should be based upon a realistic nucleon-nucleon force, since our goal is to 
obtain a unified description of the nucleon-nucleon, nucleon-nucleus and nucleus-nucleus scattering (a discussion about 
the "realism" of the interaction is found in Refs. 0||). For instance, a realistic interaction should match the empirical 
values for the volume integral and root-mean-square radius of the nucleon-nucleon interaction, Vq w —430 MeV fm 3 
and r rms « 1.5 fm, that were extrapolated from the main features of the optical potential for the nucleon-nucleus 
scattering at E nuc i eon — 10 MeV fj],46| The M3Y interaction has been derived Q with basis on the G- matrix 
for two nucleons bound near the Fermi surface, and certainly is representative of realistic interactions. In table 3 are 
presented the volume integral and root-mean-square radius for several nucleon-nucleon interactions used in this work, 
including the M3Y at 10 MeV /nucleon. 

The M3Y interaction is not truly appropriate for use in the context of the nonlocal model, because it already contains 
a simulation of the exchange effects included in its knock-on term. Furthermore, according to the nonlocal model the 



energy-dependence of the local-equivalent potential should be related only to the finite range of the Pauli nonlocality, 
but the knock-on exchange term in the M3Y interaction is also energy-dependent. Therefore, the use of the M3Y in 
the nonlocal model would imply a double counting of the energy-dependence that arises from exchange effects. In 
Section 4, we have demonstrated that the LE potential is identical with the double-folding potential for energies near 
the barrier, which are in a region around 10 MeV/nucleon. In this same energy range, the folding potential with the 
M3Y interaction have provided a very good description of elastic scattering data for several heavy-ion systems [Q. 
Thus, we believe that an appropriate nucleon-nucleon interaction for the nonlocal model could be the M3Y "frozen" 
at 10 MeV/nucleon ]l3] ], i.e. considering the parameters of the Reid and Paris versions as energy-independent values. 
Fig. 11 (top) shows a comparison between experimental and theoretical heavy-ion reduced potentials, in which the 
"frozen" M3Y-Reid was considered for the nucleon-nucleon interaction. We emphasize that no adjustable parameter 
has been used in these calculations, but even so a good agreement between data and theoretical prediction has been 
obtained. The "frozen" M3Y-Paris provides similar results. 

With the aim of investigating how much information about the effective nucleon-nucleon interaction can be ex- 
tracted from our heavy-ion potential systematics, we have considered other possible functional forms for this effective 
interaction. Besides the Yukawa function (Eq. 25), we have also used the Gaussian (Eq. 42) and the exponential (Eq. 
43), which reduce to the tri-dimensional delta function in the limit a — > 0, 

G 'W = Vo (^7^ ' (42) 



Mr)=V ^. (43) 

The fits obtained with all these functions are of similar quality and comparable with that for the M3Y interaction 
(Fig. 11 - top). The resulting best fit widths (a), volume integrals and corresponding root-mean-square radii are 
found in table 3. All the Vb and r rms values, including those of the M3Y, are quite similar. Also the "experimentally" 
extracted intensity of the nucleon-nucleon interaction in the region 1 < r < 3 fm seems to be rather independent of 
the model assumed for this interaction (see Fig. 12). 

In Section 5, we have demonstrated that the major part of the "finite-range" of the heavy-ion nuclear potential is 
related only to the spatial extent of the nuclei. In fact, even considering a zero-range for the interaction Ujviv in Eq. 
8, the shape of the heavy-ion potential could be well described just by folding the matter densities of the two nuclei. 
One would ask whether the finite-range shape of the effective nucleon-nucleon interaction can be derived in a similar 
way. Thus, we have considered a folding-type effective nucleon-nucleon interaction built from: 



v NN {r) ~ Vf(r) = I p m (ri) p m (r 2 ) V S(R - rj + f 2 ) df[ dfh = IlHYSL I n p m (n) 



n 



r 2 p m (r 2 ) dr 2 

\r-ri 



(44) 



where Vb = —456 MeV /to 3 as determined by the heavy-ion potential analysis, and p m is the matter density of the 
nucleon. Based on the intrinsic charge distribution of the proton in free space, which has been determined by electron 
scattering experiments, we have assumed an exponential shape for the matter density of the nucleon 

Pm(r) = p e~ r ' a ™ . (45) 

Of course, po and a m are connected by the normalization condition, Eq. 2. The integration of Eq. 44 results in 

«/(r) = V tt a 3 m pi e- r > a ™ ( 1 + — + £-) . (46) 

With this finite-range folding-type effective nucleon-nucleon interaction, a good fit of the reduced heavy-ion potential 
strengths is obtained (see Fig. 11 - bottom), with realistic volume integral and root-mean-square radius (see table 3). 
The folding-type interaction is quite similar to both versions of the M3Y interaction in the surface region (see Fig. 

The folding-type interaction in the context of the nonlocal model provides a very interesting unification between 
the descriptions of the nucleus- nucleus, nucleon-nucleus and effective nucleon-nucleon interactions. This can be 
appreciated through the comparison between Eqs. 36 and 38, with the subtle detail that Vf (in Eq. 36) and Vf (in 
Eq. 38) can both be calculated by folding the matter densities in the zero-range approach, and with the same Vb 
value. Therefore, the interaction between two nuclei (or nucleons) can be obtained from 



V LE (R) = / Pi(ri) P2 {r 2 ) V 5(R-fi + f 2 ) e^ 2 ^ dr[ df 2 



(47) 



where Vb = —456 MeV fm 3 , pi are the matter densities, and v is the relative speed between the nuclei (or nucleons). 
An alternative way to calculate the heavy- ion interaction is with the Eq. 37 (and 38), but in this case the nucleon 
distributions must be used (in Eq. 37) instead of the matter densities. All these findings seems to be quite consistent. 
However, the best fit value obtained for the diffuseness (a m — 0.30 fm) of the matter density of the nucleon inside 
the nucleus is considerable greater than that (a c ^ p = 0.235 fm) found for the charge distribution of the proton in 
free space. This finding is consistent with the swelling of the nucleon observed in the EMC effect p9| , but should be 
contrasted with the opposite picture of a smaller nucleon inside the nucleus as advanced within the concept of color 
transparency |5C| ]. 

Finally, we mention that, if the energy-dependence of the Pauli nonlocality is not taken into account and the 
experimental potential strengths are associated directly with the folding potential, our calculations indicate that the 
corresponding effective nucleon-nucleon interaction should have the following unrealistic values: Vq ~ —270 M eV fm 3 
and r rms ss 1.9 fm. 



7. Conclusion 

The experimental potential strengths considered in the present systematics have been obtained at the corresponding 
sensitivity radii, a region where the nuclear potential is determined from the data analyses with the smallest degree 
of ambiguity. The Fermi distribution was assumed to represent the nuclear densities, with parameters consistent with 
an extensive amount of theoretical (DHB calculations) and experimental (electron scattering experiments) results. 
The potential data set is well described in the context of the nonlocal model, by the double-folding potential in the 
zero-range as well as in the finite-range approaches. The dispersion of the potential data around the theoretical 
prediction is 25%, which is compatible with the expected effects arising from the variation of the densities due to the 
structure of the nuclei. If the nonlocal interaction is assumed, the heavy-ion potential data set seems to determine a 
few characteristics of the effective nucleon-nucleon interaction, such as volume integral and root-mean-square radius, 
in a model-independent way. 

The description of the bare potential presented in this work is based only on two fundamental ideas: the folding 
model and the Pauli nonlocality. We have avoided as much as possible the use of adjustable parameters, and in 
the case of the "frozen" M3Y interaction no adjustable parameters were necessary to fit the experimental potential 
strengths. Nowadays, the other important part of the heavy-ion interaction, the polarization potential, is commonly 
treated within a phenomenological approach, with several adjustable parameters which usually are energy-dependent 
and vary significantly from system to system. The association of the nonlocal bare potential presented in this work 
with a more fundamental treatment of the polarization should be the next step toward a global description of the 
nucleus-nucleus interaction. 

This work was partially supported by Financiadora de Estudos e Projetos (FINEP), Fundagao de Amparo a Pesquisa 
do Estado de Sao Paulo (FAPESP), and Conselho Nacional de Desenvolvimento Cienti'fico e Tecnologico (CNPq). 
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Table 1: The table presents the systems, sub-barrier bombarding energies, and corresponding references, that have 
been included in the nuclear potential systematics. 



System 


rjLab (MeV ) 


Reference 


U + INI 


OO, OO.O, OO, OO.O, O/, OO 




m\ 


KJ -f- INI 


00, OO.O, oD, 01 , 00 


I 


3C 




IbQ + b .,t>4 Ni 


34, 35, 36 




31 


i 


16Q + 88 gr 


43, 44, 45 


I 


3z 


i 


IbQ + 9U Zr 


46, 47, 48 


I 


3S 


i 


IbQ + 9^ Zr 


45, 46, 47, 48 


I 


S, 


32 


ib O + " 2 Mo 


48, 48.5, 49 


I 


32 


ib O + i2U Sn 


53, 54, 55 


I 


S 




ib + u» Ba 


54, 55, 56, 57 


I 


s 




lb + !W8p b 


74, 75, 76, 77, 78 








18 + b« M 


35.1, 35.5, 37.1, 38 


I 


Tc 




18Q + 60 M 


34.5, 35.5, 37.1, 38 


I 


3; 





Table 2: The same of table 1, but for intermediate energies. 



System 


E Lab (MeV) 


Reference 


p + 4U Ca, 2US Pb 


30.3 


I 


34 




d + 4U Ca, m Pb 


52 


I 


34 




4 He + 4U Ca, m Pb 


104 


I 


34 




b Li + i2 C, 2 «Si 


210, 318 


I 


35 


36 


y Li + 4U Ca, b «Ni, yo Zr, ™Pb 


210 


I 


37 




v Li + iy C, y «Si 


350 


I 


38 






300, 360, 1016, 1440, 2400 


I 


39 




1H C + 2U8p b 


1440 




11 




i3 C + 208p b 


390 


I 


10 




10 + ib 


250, 350, 480, 704, 1120 


I 


13 


14| 


1« + 12 Cj 28 Sij 4U Ca5 9U Zrj 2U8p b 


1504 


I 


12 




4U Ar + bU Nij laog 2U8p b 


1760 


I 


15 





Table 3: The width, volume integral and root-mean-square radius for several effective nucleon-nucleon interactions 
considered in this work. 



Interaction 


cr or a m (/to) 


V (MeV /to 3 ) 




M3Y-Rcid 




- 408 


1.62 


M3Y-Paris 




- 447 


1.60 


Yukawa 


0.58 


- 439 


1.42 


Gaussian 


0.90 


- 448 


1.56 


Exponential 


0.43 


- 443 


1.49 


Folding-type 


0.30 


- 456 


1.47 



FIG. 1. Nucleon density for the 56 Fe nucleus represented through Dirac-Hartree-Bogoliubov calculations (DHB) and a 
two-parameter Fermi distribution (2pF), with a = 0.5 fm and Ro = 4.17 fm. The small difference between the 2pF distribution 
and the function p C (Eqs. 12, 13 and 14) is hardly seen in the figure. 



FIG. 2. Equivalent diffuseness values obtained for charge distributions extracted from electron scattering experiments and 
for theoretical densities obtained from Dirac-Hartree-Bogoliubov calculations. 



FIG. 3. The Rq parameter obtained for charge distributions extracted from electron scattering experiments and for theoretical 
densities obtained from Dirac-Hartree-Bogoliubov calculations. 



FIG. 4. Nucleon densities from Dirac-Hartree-Bogoliubov calculations (solid lines) compared with the corresponding 
two-parameter Fermi distributions (dashed lines), with a — 0.50 fm and Ro obtained through Eq. 6. 



FIG. 5. Folding potential for different sets of 2pF densities that may represent the 1 O + Ni system. The approximate 
position of the s-wave barrier radius (Rb) is indicated in the figure. 



FIG. 6. Folding potential in the zero-range approach calculated from numerical integration of Eq. 8 (solid line), for 2pF 
densities that may represent the 16 O + 58 Ni system. The dashed line represents the approximate analytical expressions, Eqs. 
15, 16 and 17, while the dotted line concerns the exact result for a = 0, Eqs. 9, 10 and 11. The approximate position of the 
s-wave barrier radius (Rb) and of the distance R = Ri + R2 (s = 0) are indicated in the figure. 



FIG. 7. Normalized folding potential Vf/(Vo7Z) in the zero-range approach (Eq. 8) as a function of the distance 
s = R — (Ri + R2), for several sets of 2pF distributions (with a — 0.50 fm) that may represent the systems indicated in 
the figure. The proximity universal function $ is also presented in arbitrary units. 



FIG. 8. Double-folding potentials for 2pF distributions with different diffuseness values (a) that may represent the ie O + 
58 Ni system. The potentials have been calculated in the zero-range approach (ZR) or with a finite-range (FR) Yukawa function 
for the effective nucleon-nucleon interaction. 



FIG. 9. Double-folding (Vf) and ^-dependent local-equivalent (Vle) potentials for the a + 58 Ni system at Ehab = 139 MeV. 
The solid line represents the approximate expression, Eq. 34, for the local-equivalent potential. 



FIG. 10. Experimental and theoretical reduced potentials in the context of the zero-range approach, with (top) or without 
(bottom) considering in the calculations the energy-dependence of the local-equivalent potential (Eq. 34) that arises from the 
Pauli nonlocality. 



FIG. 11. Comparison between experimental and theoretical reduced potentials in the context of the finite-range approach, 
with a M3Y-Reid (top) or folding-type (bottom) effective nucleon-nucleon interaction. 



FIG. 12. The complete set ol effective nucleon-nucleon interactions considered in this work. 
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